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ABSTRACT 


Stanford  Research  Institute,  under  the  sponsorship  of  ARPA 
and  Office  of  Naval  Research,  is  studying  two  dif luoroamino  com¬ 
pounds  in  the  following  three  situations. 

A.  The  detonation  characteristics,  including  the 
shock  sensitivity,  of  the  liquid  phase 

B.  The  adiabatic  self -heating  of  the  liquid  phase 

C.  The  mechanism  and  kinetics  of  thermal  decomposition 
of  the  liquid  phase. 

These  approaches  were  chosen  for  their  susceptibility  to 
analysis  in  terms  of  basic  physical  and  chemical  parameters  and 
because  they  are  interrelated.  Each  study  is  being  carried  out 
on  the  isomers,  1 , 2-bis (dif luoroamino)propane  and  2 , 2-bis(dif luoro- 
amino)propane ,  in  order  to  assess  the  importance  of  structural 
variation . 

A.  Detonation  Characteristics.  The  objectives  of  this  phase  of  the 
work  are  (l)  to  determine  quantitatively  some  of  the  detonation  charac¬ 
teristics  of  the  two  compounds,  1,2- and  2 , 2-bis (dif luoroamino) propane , 
including  the  kinds  of  shock  that  will  initiate  the  detonation,  and  (2) 
to  study  the  relationship  of  shock  sensitivity  to  the  other  detonation 
characteristics,  to  the  physical  properties,  and  to  the  chemical  kinetics 
of  thermal  decomposition.  Experimental  work  during  the  past  year  has 
included  measurements  of  detonabil ity ,  detonation  velocity,  and  failure 
diameter  of  1,2  and  2,2  DP.  Both  cyclindricai  and  conical  confinements 
were  used  for  the  experiments  with  the  1,2  DP  and  conical  confinements 
for  2,2  DP.  The  conical  confinement  consisted  of  a  combination  of 
straight,  thin-walled  copper  tube  (long  enough  to  obtain  stable  detona¬ 
tion)  with  tapering  to  about  1-cm  diameter  over  the  rest  of  the  length. 

The  experiments  establish  that  both  1,2  and  2,2  DP  detonate  at 
velocities  of  5-7  mm/usec  and  have  failure  d/ 'raters  of  less  than  10  mm. 
Variations  of  detonation  velocity  from  shot  to  shot  and  within  one  shot 
are  greater  than  can  be  accounted  for  by  the  uncertainties  in  the  measuring 
techniques.  This  variation  is  as  yet  unexplained. 

The  ideal  Chapman-Jouguet  values  were  calculated,  using  the  B.K.W. 
equation  of  state  and  the  Ruby  Code  method,  for  the  reported  density 
value  of  1.256  g/cc;  a  later  measurement  in  our  laboratory  gave  1.17 
g/cc.  The  C-J  velocity  was  6.55  mm/fosec;  wie  pressure  was  192  Kbar;  the 
temperature  was  2828°K. 
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B.  Adiabatic  Self-Heating.  If  an  explosive  is  placed  in  an  enviror 
ment,  the  temperature  of  which  is  continuously  changed  so  as  to  follow 
precisely  the  rising  temperature  of  the  decomposing  sample,  then  a  modi¬ 
fication  of  the  Arrhenius  rate  equation  for  a  first  order  reaction  may 

be  applied  in  its  logarithmic  form: 

ln(||)  =  ln(qZ/c)  -  E/R (l/T ) 

Hence,  given  a  temperature-time  history  from  a  single  experiment  and  a 
value  of  the  heat  capacity,  one  may  obtain  values  qZ  and  E.  For  a  fluid 
requiring  a  confining  vessel  the  mathematical  treatment,  though  more 
complicated,  is  still  tractable  in  terms  of  calculating  E  and,  with  cer¬ 
tain  assumptions,  qZ . 

Self-heating  experiments  of  this  type  have  been  conducted  on  two 
liquid  compounds:  1,2  DP  and  2,2  DP.  The  results  indicate  that  the 
thermal  decomposition  of  1,2  DP  does  not  proceed  as  required  by  the 
mathematical  model,  giving  an  unrealistic  energy  of  activation  value. 

On  the  other  hand  the  2,2  DP  appears  to  decompose  in  a  simple  fashion. 

The  energy  of  activation  calculated  from  three  runs  is  AEa  =  51.0  ±  2.5 
kcal/mole.  A  value  of  qZ  will  be  calculated  when  the  additional  required 
data  are  obtained. 

C.  Kinetics  and  Mechanism  of  Thermal  Decomposition.  The  thermal 
decomposition  of  1,2  DP  has  been  examined  in  nitrobenzene  solution  over 
the  range  176-187°.  The  initial  rate  was  found  to  be  invariant,  within 
experimental  error,  over  a  fourfold  change  in  concentration  (0.52-0.14  M) . 
After  one  quarter  to  one  half-life,  the  rate  markedly  accelerates.  The 
following  rate  laws  were  found  to  describe  the  reaction  over  much  of  the 
decomposition : 

Rate  (176«>)  =  dx/dt  =  k,(A-x)  +  80ki(A-x)x2;  k,  =  1  x  10_s  sec"1 
Rate  (187°)  =  dx/dt  =  k,(A-x)  +  64k,(A-x)x2;  k,  =  3  x  10-5  sec-1 

Kinetic  parameters  for  both  the  first-order  portion  and  autocatalytic 
portion  of  the  decomposition  were  estimated: 

k,  =  3.7  x  1014  exp (-40 .2  ±  2/RT)  sec-1  (Ea  =  kcal/mole 
k3  =  1.8  x  1014  exp ( -35 .7  ±  2/RT)  l2  mol-2  sec"1 

The  principal  reaction  product  from  long  heating  of  1,2  DP  is  a  black, 
insoluble  polymer  of  unknown  structure.  Several  volatile  intermediates 
have  been  detected  in  small  amounts.  Two  of  these  intermediates  are 
thought  to  be  acetonitrile  and  2- (N-f luoroimino)propionitrile .  This 
same  pattern  of  products  is  observed  in  the  decomposition  of  both  the 
pure  liquid  and  nitrobenzene  solutions. 

The  rate  of  decomposition  of  1,2  DP  in  nitrobenzene  is  markedly 
increased  by  small  amounts  of  2, 6 -di-t -butyl -4 -methyl phenol ,  a  free- 
radical  chain  inhibitor.  a-Methyl styrene  also  Increases  the  rate  of 
decomposition.  Water  catalyzes  the  decomposition,  but  does  not  appear 
to  be  responsible  for  the  observed  autocatalysis.  Replacement  of 
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nitrobenzene  by  hexadecane  leads  to  a  diminution  in  the  rate  of  decom¬ 
position.  Sodium  butoxide  in  butanol  destroys  1,2  DP  rapidly  at  25°; 
the  volatile  intermediates  produced  in  this  reaction  are  identical  to 
those  produced  in  the  thermal  decomposition. 

Pronounced  surface  effects  were  found  in  the  vapor-phase  decomposi¬ 
tion  of  1,2  DP  (10~3  M  at  187°).  The  reaction  is  accelerated  in  steel 
vessels  but  retarded  in  glass  vessels.  In  a  new  glass  bulb  the  rate  is 
about  one-tenth  the  rate  found  in  nitrobenzene  solution.  Products  of  the 
vapor-phase  decomposition  are  identical  with  those  found  in  the  solution 
reaction;  this  suggests  that  the  vapor-phase  decomposition  takes  place 
at  the  surface. 

We  believe  that  the  rate-controlling  step  in  the  thermal  decomposi¬ 
tion  of  1,2  DP  involves  the  loss  of  HF  from  the  secondary  carbon,  followed 
by  C-C  bond  cleavage  and  polymerization  of  the  intermediates.  Probably, 
the  high  molecular  weight  residue  is  responsible  for  the  observed  auto¬ 
catalysis  . 

The  vapor  phase  decomposition  of  2,2  DP  was  examined  at  200°  in  glass 
and  steel  bulbs  at  10“3  M  concentration.  The  decomposition  was  found  to 
be  autocatalytic .  Several  intermediate  products  were  detected  but  the 
final  product  is  a  dark,  polymeric  residue.  The  rate  of  decomposition 
is  about  one-tenth  that  of  1,2  DP.  Probably  the  reaction  occurs  on  the 
surface . 

A  brief  examination  of  the  decomposition  of  2,2  DP  in  nitrobenzene 
indicates  the  rate  here  is  also  much  slower  than  the  rate  for  1,2  DP. 
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I  INTRODUCTION 

The  word  sensitivity  implies  a  quantitative  measure  of  the  ease  of 
inducing  in  a  material  a  rapid  exothermic  chemical  decomposition  accom¬ 
panied  by  a  pressure  rise  of  sufficient  magnitude  to  be  destructive. 

Many  empirical  tests  devised  to  measure  sensitivity  fail  to  correlate, 
are  only  partially  satisfactory  as  predictive  tests,  and  give  little 
fundamental  information.  Moreover,  the  results  of  such  tests  are  not 
interpretable  in  terms  of  basic  physical  and  chemical  behavior. 

There  are  two  reasons  for  this  unsatisfactory  state  of  affairs. 

First,  many  of  the  tests  are  so  complicated  that  they  cannot  be  reduced 
to  the  basic  physics  and  chemistry  of  the  system.  From  the  results,  one 
cannot  analyze  the  physical  events  reliably  enough  to  deduce  the  chemis¬ 
try;  nor  does  one  know  the  chemistry  well  enough  to  be  able  to  interpret 
the  results  in  terms  of  the  applied  disturbance.  Second,  there  are 
several  different  kinds  of  events  the  results  of  which  can  be  termed 
explosive.  For  example,  a  given  high  energy  liquid  chemical,  depending 
upon  circumstances  and  upon  its  chemical  and  physical  characteristics, 
may  undergo  one  or  more  of  the  following  kinds  of  exothermic  reaction 
accompanied  by  a  very  fast  pressure  rise:  a  detonation  wave  in  the  liquid 
phase;  thermal  decomposition  in  the  liquid  phase  uniformly  throughout 
the  mass;  vaporization  of  the  liquid,  followed  by  a  detonation  wave  in 
the  gas  phase;  vaporization  of  the  liquid,  followed  by  thermal  decomposi¬ 
tion  in  the  gas  phase;  vaporization  of  the  liquid,  followed  by  initiation 
of  a  deflagration  wave  in  the  gas  phase,  and  then  by  ignition  and  defla¬ 
gration  in  the  liquid  phase. 

The  reason  for  the  lack  of  correlation  of  the  empirical  sensitivity 
tests  may  often  be  that  the  tests  are  concerned  with  different  events 
in  the  lists  above,  i.e.,  different  paths  of  decomposition.  A  card-gap 
test  measures  sensitivity  to  detonation  of  the  liquid  phase.  A  drop- 
weight  or  impact  test  may  reflect,  at  least  in  part,  the  ease  of  thermal 
decomposition  of  the  liquid  phase.  The  hot-wire  test  may  measure  the 
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ease  with  which  the  liquid  vaporizes  and  ignites  in  the  gas  phase. 

Thermal  stability  measurements  may  reflect  liquid  phase  decomposition 
kinetics . 

Clearly,  it  is  desirable  to  define  precisely  the  event  during  any 
sensitivity  measurement.  In  order  to  predict  with  assurance  the  kind 
of  explosive  behavior  to  be  expected  of  a  material  under  a  given  set  of 
circumstances,  one  should  know: 

1.  Thermal  stability  in  terms  of  liquid  phase  decomposition 
kinetics 

2.  Detonability  of  liquid  phase,  i.e.,  whether  it  will  support 
a  detonation  wave,  and  the  failure  diameter  for  detonation 

3.  The  necessary  conditions  for  initiating  a  detonation  wave 
in  the  liquid  phase 

4.  The  def lagrability  of  the  liquid  phase,  i.e.,  the  conditions 
of  pressure  and  geometry,  if  any,  under  which  the  liquid 
supports  a  deflagration  wave  at  its  surface. 

5.  The  necessary  conditions  for  initiating  a  deflagration  wave 
at  the  liquid  surface 

6.  The  decomposition  kinetics  of  the  vapor  phase 

7.  Detonability  of  the  vapor  phase 

8.  The  necessary  conditions  for  initiating  a  detonation  wave 
in  the  vapor  phase 

9.  The  deflagrability  of  the  vapor  phase 

10.  The  necessary  conditions  for  initiating  a  deflagration  wave 
in  the  vapor  phase. 

A  common  need  for  understanding  each  of  the  ten  decomposition  paths 
is  knowledge  of  the  decomposition  chemical  kinetics.  In  this  work  our 
aim  is  to  study  the  decomposition  kinetics  and  their  relationship  to  one 
specific  sensitivity  situation  which  can  be  relatively  clearly  defined, 
namely,  the  shock  sensitivity.  We  expect  the  information  on  decomposi¬ 
tion  kinetics  to  be  useful  in  understanding  and  predicting  sensitivity  in 
other  experimental  situations  as  well. 
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We  have  undertaken  the  study  for  several  liquid  materials  of  re¬ 
lated  chemical  structure,  hydrocarbons  containing  two  dif luoroamino  (NF2 ) 
groups.  The  compounds  being  studied  initially  are  two  isomers,  1,2- 
bis(dif luoroamino)propane  and  2 , 2-bis (dif luoroamino)propane .  These  two 
compounds  are  hereafter  referred  to  as  1,2  DP  and  2,2  DP,  respectively. 
The  program  is  divided  into  three  interrelated  parts.  They  are 

(a)  The  shock  sensitivity  to  detonation  of  the  liquid  phase 

(b)  The  adiabatic  self-heating  of  the  liquid  phase 

(c)  The  mechanism  and  kinetics  of  thermal  decompositior 

In  (a)  the  detonation  characteristics,  particularly  the  shock  sensitivity 
to  detonation,  are  measured.  The  results  will  be  analyzed  in  terms  of 
detonation  theory,  taking  into  account  the  physical  properties  of  the 
liquid  and  the  chemical  kinetics  of  the  liquid  phase  decomposition. 
Adiabatic  self-heating,  (b),  is  a  method  for  obtaining  values  for  the 
pre-exponential  factor  and  the  activation  energy  of  the  liquid  phase 
decomposition  if  the  over-all  reaction  rate  can  be  formally  treated  as 
a  simple  first-order  rate  process.  The  simple  form  of  the  rate  equation 
thus  obtained  makes  it  attractive  for  use  in  theories  of  detonation, 
critical  size  for  ignition,  thermal  explosion,  and  impact  sensitivity. 
Part  (c),  in  addition  to  providing  rate  data  for  various  theories  of 
sensitivity,  will  provide  information  assisting  in  an  understanding  of 
the  relationship  between  chemical  structure  and  sensitivity. 
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II  SHOCK  SENSITIVITY 
(Marjorie  W.  Evans  and  A.  B.  Amster) 


A .  Summary 

The  objectives  of  this  phase  of  the  work  are  to  determine  quanti¬ 
tatively  some  of  the  detonation  characteristics  of  1,2  and  2,2  DP,  in¬ 
cluding  the  shock  sensitivity  to  detonation,  and  to  study  the  relationship 
of  shock  sensitivity  to  the  other  detonation  characteristics,  to  the 
physical  properties,  and  to  the  chemical  kinetics  of  thermal  decomposit ior 

Experimental  work  during  the  past  year  has  included  measurements  of 
detonabil ity ,  detonation  velocity,  and  failure  diameter  of  1,2  DP  and 
2,2  DP.  Theoretical  work  has  included  a  calculation  of  the  ideal  Chapman- 
Jouguet  detonation  velocity,  pressure  and  temperature  of  1,2  DP,  using 
a  method  which  has  been  reasonably  successful  in  calculating  these 
quantities  in  CHNOF  systems.  In  addition,  the  equations  describing  detona 
tion  characteristics  of  homogeneous  condensed  materials  in  terms  of  their 
chemical  and  physical  properties1  have  been  programmed  and  solutions  are 
being  obtained  for  a  range  of  parameters. 

B .  Background 

Shock  sensitivity  is  determined  bj  measuring  the  characteristics 
of  a  minimum  initiating  shock  which  is  moving  into  the  explosive  from 
an  inert  attenuating  material .  The  sense  in  which  "minimum"  is  used 
requires  definition.  The  shock  variables  of  pressure,  density,  and 
material  velocity  in  the  general  case  are  space-  and  time-dependent. 

For  convenience  in  treatment  it  is  helpful  to  simplify  the  shock  or  to 
limit  the  space-time  region  which  is  said  to  influence  the  initiation  of 
detonation,  or  to  do  both.  As  an  example  of  the  first,  one  can  specify 
plane  shocks  of  square  profile  and  define  shock  sensitivity  In  terms  of 


‘Marjorie  W.  Evans,  J.  Chem.  Phys.  36,  193  (1962) 
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magnitude  and  duration  of  initiating  shock.  As  an  example  of  the  second, 
one  can  make  the  assumption  that  in  cylindrical  charges  the  controlling 
events  occur  on  the  axis,  where  the  shock  front  is  perpendicular  to  the 
axis  and  the  shock  pressure  is  a  maximum;  then  one  can  substitute  for 
the  planar  shock  the  more  easily  created  curved  shock,  and  measure  the 
shock  characteristics  on  the  axis.  With  the  approximation  that  the  shock 
is  square,  the  shock  pressure  on  the  axis  at  the  initiating  boundary  of 
the  explosive  rises  rapidly  to  a  value  which  remains  constant  for  a  dura¬ 
tion  t,  at  which  time  a  rarefaction  wave  from  side  or  rear  abruptly 
reduces  the  pressure  to  the  ambient  value.  For  this  latter  situation 
the  minimum  shock  sensitivity  is  conveniently  defined,  for  shocks  of 
given  duration,  in  terms  of  the  minimum  pressure,  p,  required  for  initia¬ 
tion.  Given  the  correctness  of  the  assumption  and  the  appropriateness  of 
the  approximation,  measurements  of  the  shock  sensitivity  by  the  two 
systems,  plane  square  shock  and  curved  shock,  will  be  in  accord.  It 
seems  clear  that  shock  sensitivity  must  be  defined  as  a  function  of  the 
two  variables  magnitude  and  duration .  When  the  duration  of  a  shock  is 
limited  by  the  diameter  of  the  receiver  rather  than  the  physical  dimension 
of  the  booster  explosive,  it  follows  that  the  minimum  magnitude  of  shock 
for  initiation  will  depend,  for  a  cylindrical  charge,  on  its  diameter. 
Thus,  if  shock  sensitivity  is  in  such  a  case  defined  as  magnitude  of 
shock,  one  expects  it  to  vary  with  diameter. 

Relations  connecting  the  detonation  characteristics  (including  shock 
sensitivity)  and  the  physical  and  chemical  characteristics  have  been 
worked  out  for  laminar  detonation  of  homogeneous  condensed  materials 
having  thermal  decomposition  rates  expressible  as  a  single  first-order 
rate  process.1  The  situation  in  that  case  is  simplified  by  two  circum¬ 
stances.  First,  the  form  of  the  rate  law  dictates,  in  a  shock  of  given 
magnitude  and  duration,  that  chemical  reaction  occurs  either  essentially 
completely  or  essentially  not  at  all.  Second,  the  mechanism  of  initiation 
of  chemical  reaction  at  the  front  of  a  formed  detonation  wave  is  the  same 
as  the  mechanism  of  initiation  of  chemical  reaction  in  a  shock,  i.e., 
each  occurs  by  shock  compression  to  a  temperature  sufficiently  high  that 
reaction  takes  place  before  a  rarefaction  wave  intervenes.  The  simplicit> 
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of  the  system  makes  possible  the  statement  of  a  set  of  algebraic  equations 
which  not  only  yield  the  dependence  of  shock  sensitivity  on  diameter,  but 
also  relate  this  in  an  internally  consistent  model  to  failure  diameter, 
induction  time,  and  the  relationship  between  diameter  and  steady  detona¬ 
tion  velocity. 

Before  measuring  the  shock  sensitivity  of  a  given  material,  one  must 
establish  its  detonabil ity ,  which  is  the  ability  of  a  material  to  support 
a  steady  detonation  wave.  It  is  a  function  of  diameter,  the  minimum 
diameter  for  which  a  material  is  detonable  being  the  failure  diameter . 
Experimental  work  during  this  year  has  been  directed  toward  studying 
the  detonability  of  1,2  and  2,2  DP.  It  is  necessary  to  establish  detona¬ 
bil  ity  with  sufficiently  long  charges  that  unequivocal  demonstration  is 
made  of  either  the  formation  of  a  steady  detonation  wave,  or  the  failure 
of  a  transient  reaction  wave. 

C .  Experimental 


Three  methods  were  used  for  determining  detonability  and  detonation 
velocity:  ionization  probes,2'3’4  pressure  probes,5  and  a  continuous 

wire  sensor6  (the  latter  may  be  referred  to  as  CRISP — continuous  recording 
of  instantaneous  shock  properties).  In  one  method  of  using  the  continu¬ 
ous  wire  sensor,  a  fine  high  resistance  wire  is  placed  in  the  cylindrical 
charge,  parallel  to  the  long  axis.  The  charge  is  enclosed  in  a  metal 
case.  A  voltage  from  a  constant  current  source  is  applied  at  the  end  of 
the  wire  opposite  the  end  at  which  the  detonation  is  initiated.  The 
ionization  in  the  moving  detonation  wave  completes  the  circuit  from  the 
wire  to  the  case.  Since  the  voltage  drop  along  the  wire  is  directly 
proportional  to  the  resistance  of  the  wire  remaining  in  the  circuit,  it 
is  also  proportional  to  the  length  of  undetonated  charge.  Thus  from  a 
measure  of  the  voltage  drop  and  its  change  with  time  one  can  calculate 
the  detonation  velocity  as  a  function  both  of  time  and  of  position  of 

2 Frank  C.  Gibson,  Rev.  Sci.  Instr.  25^  226  (1954) 

3H.  T.  Knight  and  R.  E.  Duff,  Rev.  Sci.  Instr.  26,  257  (1955) 

4A.  W.  Campbell,  M.  E.  Malin,  T.  J.  Boyd,  Jr.,  and  J.  A.  Hull,  Rev.  Sci. 

Instru.  27,  567  (1956). 

*A .  B.  Amster  and  R.  L.  Beauregard,  Rev.  Sci.  Instr.  30,  942  (1959) 

eA.  B.  Amster,  P.  A.  Kendall,  L.  J.  Veillette,  and  Burton  Harrell,  Rev. 

Sci.  Instr.  31,  188  (i960) 
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the  detonation  front.  A  typical  record  of  voltage  drop  versus  time,  as 
displayed  on  a  cathode  ray  oscilloscope,  is  shown  in  Fig.  1.  The  oscil¬ 
loscope  is  triggered  at  time  A,  and  the  wave  reaches  the  wire  at  time  B. 
The  line  CD  records  the  changing  voltage  drop  as  the  wave  moves  parallel 
to  the  wire.  If  the  detonation  velocity  is  constant,  CD  is  straight. 
When  detonation  ceases,  the  voltage  drops  to  E,  Lines  I  and  II  are  in¬ 
dependently  recorded  to  provide  time  and  voltage  calibrations. 


FIG.  1  TYPICAL  OSCILLOSCOPE  RECORD  OF  CONTINUOUS  WIRE 
DETONATION  VELOCITY  MEASUREMENT 


The  1,2  and  2,2  DP  were  prepared  immediately  before  the  detonation 
tests  by  distillation  from  the  solvent  (see  the  Appendix).  Chromatographic 
analyses  showed  that  the  material  was  99.4$  pure  or  better.  Three  tests 
each  were  performed  on  1,2  and  2,2  DP. 

Four  experimental  arrangements  were  used,  as  shown  in  Figs.  2,  3, 

4,  and  5.  The  cylindrical  arrangement  of  Fig.  2  corresponds  to  the  card- 
gap  test  recommended  by  the  JANAF  Panel  on  Liquid  Propellant  Test  Methods 
except  for  the  presence  of  ionization  probes  and  the  absence  of  a  witness 
plate.  The  cylindrical  arrangement  of  Fig.  3  incorporated  the  simultan¬ 
eous  use  of  ionization  probes,  pressure  probes,  and  the  continuous  wire 
sensor.  In  the  arrangements  of  Figs.  4  and  5  the  cylindrical  section 
provides  a  length  in  which  detonation  can  achieve  a  steady  velocity; 
failure  diameter,  if  it  is  greater  than  the  minimum  diameter  of  the  cone, 
will  be  determined  by  the  position  of  the  failure  in  the  conic  section. 
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^  l  INCH  SCHEDULE  40 
BLACK  STEEL  PIPE 


FIG.  2  DETONATION  VELOCITY  APPARATUS  (I) 
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CROSS  SECTION  OF  FRONT  VIEW 

A. 40||  .lid 


FIG.  3  DETONATION  VELOCITY  APPARATUS  (II) 
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FIG.  4  DETONATION  VELOCITY  APPARATUS  (III) 
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FIG.  5  DETONATION  VELOCITY  APPARATUS  (Ill-Modified) 


One  test  with  the  1,2  DP  was  carried  out  in  each  of  three  experimental 
arrangements.  All  three  tests  with  the  2,2  DP  were  in  the  arrangement  of 
Fig.  5.  The  observed  failure  diameters  and  the  detonation  wave  velocities 
as  a  function  of  distance  from  the  attenuator-acceptor  Interface  are 
given  in  Table  I.  Both  distances  and  failure  diameters  are  calculated 
from  and  must  reflect  the  precision  of  the  velocities.  The  experiments 
establish  that  1,2  and  2,2  DP  detonate  at  velocities  of  5-7  mm/nsec,  and 
have  failure  diameters  of  less  than  10  mm.  The  variations  of  velocity 
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from  shot  to  shot  and  within  one  shot  are  greater  than  can  be  accounted 
for  by  the  uncertainties  in  the  measuring  techniques.  We  cannot  for  the 
present  explain  the  discrepancy. 


O .  Theoretical 

The  ideal  Chapman-Jouguet  detonation  characteristics  were  calculated 
for  1,2  DP  using  a  method  which  has  been  reasonably  successful  in  calcu¬ 
lating  detonation  characteristics  of  CHNOF  systems.*  It  is  based  on  the 
Kistiakowsky-Wilson  equation  of  state,  as  modified  by  Cowan  and  Fickett7 
and  generalized  by  Mader.8  The  Brinkley  method  is  used  for  computing  the 
equilibrium  composition,  and  the  generalized  technique  is  referred  to 
as  the  BKW  calculation.  The  calculation  was  made  using  the  Ruby  Code.9 
It  permits  computing  the  equilibrium  composition  of  the  explosion  products 
the  detonation  Hugoniot ,  the  values  of  the  hydrodynamic  and  thermodynamic 
variables  at  the  C-J  point,  and  the  isentrope  of  the  reaction  products. 
Table  II,  which  gives  input  data  and  results,  reports  the  calculated 
equilibrium  composition  of  products  and  the  values  of  the  C-J  pressure, 
temperature,  and  detonation  velocity. 

The  covolumes  of  the  products  and  the  constants  in  the  equation  of 
state  were  those  suggested  by  Mader.8  The  required  reference-state 
thermodynamic  functions  and  enthalpies  of  formation  were  taken  from  the 
JANAF  Interim  Chemical  Tables.10  The  calculation  was  performed  neglecting 
CH4  in  the  product  gas.  It  is  not  expected  that  the  values  of  the  C-J 


*The  calculations  were  made  by  personnel  of  the  Chemistry  Division  at 
Lawrence  Radiation  Laboratory,  Livermore,  California.  We  are  particu¬ 
larly  indebted  to  Dr.  Edward  Lee  for  his  help. 

7R.  D.  Cowan  and  W.  Fickett,  J.  Chem.  Phys .  24,  932  (1956) 

8C .  R.  Mader,  Los  Alamos  Scientific  Laboratory  Report  LA-2613,  Oct. 

9,  1961,  "Detonation  Performance  Calculations  Using  the  Kistiakowsky- 
Wilson  Equation  of  State." 

9H.  B.  Levine  and  Robert  E.  Sharpies,  Lawrence  Radiation  Laboratory 
Report  UCRL-6815,  March  20,  1962,  "Operator's  Manual  for  Ruby." 

10 JANAF  Interim  Chemical  Tables,  The  Dow  Chemical  Company,  Midland, 
Michigan . 
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Table  II 

INPUT  DATA  AND  RESULTS  OF  DETONATION  PERFORMANCE 
CALCULATION  FOR  1,2  DP 


Formula 

Mol  Wt 

Density  (g/ccj 

(kcal/mole  HE) 

Tc  (°K) 

i 

c3h6f4n2 

146.09 

1.256* 

-40.1** 

298 

C-J  Gas  Volume  (cc/mole 

product  gas) 

17.7 

Moles  Product  Gas/gram 

reactant 

0,0253 

Moles  Solid  Carbon/gram 

reactant 

0.0156 

Gamma 

1,808 

No.  Moles  Product  Gas/gram  reactant 

h2 

0.00029 

NF 

0.0 

nh3 

0.01086 

cf4 

0.00489 

n2 

0.00142 

HF 

0.00784 

Calculated  C-J  values 

C-J  Pressure  (kbar)  192 

C-J  Temperature  (°K)  2828 

Detonation  Velocity  (mm/usec)  6.55 


*First  Quarterly  Progress  Report  on  Advanced  Test  Methods  for 
Determining  Operational  Characteristics  of  Propellants, 
Document  No  61-21-16  Contract  AF  04(611 )-7413 ,  Air  Reduction 
Company 

♦♦Reported  at  the  July  1962  JANAF  Thermochemical  Panel. 


pressure,  temperature,  and  detonation  velocity  will  be  greatly  affected 
if  it  is  taken  into  account.  The  value  of  gamma  is  low,  possibly  as  a 
result  of  the  relatively  low  value  for  moles  of  product  gas  per  gram  of 
reactant.  The  C-J  pressure  and  temperature  are  in  the  range  expected. 
Since  these  calculations  were  made  a  measurement  of  density  in  our  labora¬ 
tories  gave  a  value  of  1.18  g/cc  at  25°C  rather  than  1.256.  Detonation 
calculations  based  on  this  lower  value  of  density  can  be  expected  to  give 
lower  values  of  detonation  velocity  (about  6  mm/pisec)  and  of  pressure. 
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The  equations1  describing  characteristics  of  two-dimensional  laminar 
detonation  waves  in  homogeneous  condensed  materials  have  been  programmed, 
and  their  solutions  are  being  obtained  for  a  range  of  parameters  which 
include  those  expected  to  be  appropriate  for  the  dif luoroamino  compounds. 

E.  Future  Plans 


The  detonation  velocities  of  1,2  and  2,2  DP  and  of  nitromethane  will 
be  studied  further  in  an  attempt  to  establish  the  significance  of  the  two 
levels  of  detonation  velocity  reported  in  Table  I.  Observational  tech¬ 
niques  will  include  the  use  of  a  streak  camera  in  addition  to  those  already 
employed.  Nitromethane  will  be  included  in  the  program  as  a  check  com¬ 
pound,  since  its  detonation  properties  are  reasonably  well  established.11" 
Measurements  of  shock  sensitivity  will  be  made.  The  detonation  velocities 
and  shock  sensitivities  of  IBA ,  1,3-  and  1 , 1-bis (dif luoroamino) propanes 
will  be  studied  as  they  become  available  and  can  be  worked  into  the 
schedule . 

The  computations  based  on  the  theoretical  equations  will  be  carried 
out  over  a  range  of  parameters  believed  to  be  suitable.  As  experimental 
data  on  physical  and  chemical  properties  and  on  detonability  and  sensi¬ 
tivity  become  available,  the  theoretically  predicted  detonation  charac¬ 
teristics  will  be  compared  with  observations,  with  a  view  to  establishing 
the  usefulness  of  the  theory  and  to  improving  it. 


1 *A .  W.  Campbell,  M.  E.  Malin,  and  T.  E.  Holland, J.  Appl .  Phys.  27, 

963  (1956). 

12A.  W.  Campbell,  W.  C.  Davis,  and  J.  R.  Travis,  Phys.  Fluids  £,  498  (1961 
13 J.  R.  Travis,  A.  W.  Campbell,  W.  C.  Davis,  and  J.  B.  Ramsay,  "Shock 
Initiation  of  Explosives  III.  Liquid  Explosives,"  Les  Ondes  de 
Detonation,  Coll.  Inter.  Centre  Nat.  Rech .  Scien.,  1962,  p.  45 
14G.  Nahmani  and  Y.  Manheimer,  J.  Chem.  Phys.  24,  1074  (1956). 
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III  ADIABATIC  SELF-HEATING 
(A.  B.  Amster) 


A .  Summary 

The  adiabatic  self-heating  method15  for  determining  activation 
energy  and  pre-exponential  factor  of  condensed  phase  thermal  decomposi¬ 
tion  was  applied  to  1,2  and  2,2  DP.  Results  with  the  1,2  DP  were  in¬ 
consistent,  and  yielded  unreasonably  high  values  of  activation  energy. 
Results  with  the  2,2  DP  were  quite  reproducible  and  gave  a  value  of  51.0 
±2.5  kcal/mole  for  the  activation  energy. 


B .  Background 


Consider  a  finite  quantity  of  material  at  temperature  T,  in  an 
environment  of  some  other  temperature  T1 .  If  the  material  decomposes 
according  to  a  simple  zero-order  process,  its  temperature  is  given  by 
the  following  equation: 

XV2T  +  pc(|£)  =  Ze"Ea//RT  (3-1) 

and  further,  if  the  mechanism  ?s  first  order,  Z  =  pQv  and 

-  aV2T  +  pc(|£)  =  PQVe  Ea^RT  (3-2) 


where 

t  =  time  (sec) 

X  =  thermal  conductivity  (cal  deg" 1  cm"1 see"1 ) 
p  =  density  (g  cm"3 ) 
c  =  specific  heat  (cal-1  g*1  deg*1) 

T  =  temperature  (°K) 
q  =  heat  of  decomposition  (cal  g"1) 

V  =  frequency  factor  (sec"1) 

Ea  =  activation  energy  (cal  IT1) 

R  =  gas  constant  (1.987  cal  M"1  deg-1) 

Z  =  heat  generation  coefficient  (cal  cm"3  sec"1) 


ISD.  Gross  and  A.  B.  Amster,  Eighth  Symposium  (international)  on  Combust ic 
Williams  and  Wilkins  Company,  Baltimore,  1962,  pp.  728-734. 
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In  the  particular  case  in  which  the  environmental  temperature  is  con¬ 
tinuously  changed  so  as  to  follow  precisely  the  rising  temperature  of 
the  decomposing  sample,  V2T  =  0  and 


or 


qve 


-Ea/RT 


1 


dT 

™dt 


Ea 

RT 


(3-3) 

(3-4) 


Hence,  given  a  temperature-time  history  from  a  single  experiment  and 
a  value  of  the  heat  capacity,  one  may  obtain  values  for  qv  and  Efl .  For 
non-first-order  reactions,  the  evaluation  of  qv  is  slightly  more  involved 
but  can  still  be  accomplished;  the  interpretation  of  Ea  is  unambiguous. 
The  case  of  simultaneous  reactions  is  considered  in  section  D  below. 


When  the  sample  material  is  a  liquid  the  experimental  arrangement 
is  somewhat  more  complicated  than  with  solids,  and  the  above  analysis 
must  be  modified.  To  prevent  evaporation,  the  sample  is  confined  within 
a  metal  vessel  capable  of  withstanding  the  vapor  pressure  of  the  liquid 
at  elevated  temperatures.  The  confining  vessel  is  now  heated  partially 
by  the  heat  generated  by  the  liquid  decomposition  and  partially  by  the 
furnace,  and  the  analysis  becomes: 


Let 

mL  =  mass  of  liquid  (g) 

mM  =  mass  of  fraction  of  container  heated  by  liquid  (q) 
The  heat  generation  rate  of  the  liquid,  dQ/dt  cal/sec,  is 


dQ 

dt 


mLqVe 


-Ea/RT 


(3-5) 


If  there  is  no  temperature  gradient  within  the  liquid  or  the  containing 
vessel  or  between  them,  the  rate  of  increase  in  heat  content  for  the 
liquid  and  that  part  of  the  container  heated  by  the  liquid  is: 


un  .  v  ui 

—  =  (m||cM  ♦  «»Lcl)  dt 


(3-6) 
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Equating  Eqs.  (3-5)  and  (3-6)  and  taking  natural  logarithms  gives 

dT  „  mLQV  Ea 

In—  =  In  -  -  — 

dt  mMcM  +  mLcL  RT 

Thus,  a  plot  of  ln(dT/dt)  vs.  l/T  should  give  a  straight  line  with  slope 
-E/R  and  with  intercept 


1  n 


mLqV 

mMcM  +  mLcL 


(3-8) 


where  c^,  c^,  and  are  known  or  can  be  evaluated  independently.  An 
evaluation  of  therefore  permits  determination  of  qv . 


An  approximate  value  for  can  be  obtained  as  follows:  Assume  that 
the  boundary  between  the  regions  in  the  metal  vessel  heated  by  the  liquid 
and  the  oven  lies  halfway  between  the  outer  and  inner  surfaces  of  the 
approximately  spherical  vessel .  If  is  the  radius  of  the  cavity  and 
r2  the  outer  radius  of  the  vessel,  then  the  mass  of  metal  heated  by  the 
liquid  is 

-  mM  (3-9) 

where  Pm  is  the  density  of  the  metal  vessel .  This  value  can  then  be  sub¬ 
stituted  in  expression  (3-8)  and  qv  calculated,  since  the  other  parameters 
are  known.  From  experiments  with  bombs  of  differing  thicknesses  the 
reliability  of  this  approximation  can  be  tested.  If  the  model  is  suitable, 
extrapolation  of  these  intercepts  to  zero  thickness  will  give 

qv 

ln  mMCM  +  m"LcL  =  ^  -  for  ^  =  0  (3-10) 


C  Experimental 

1 .  Experimental  Details 

The  equipment  consisted  of  (a)  pressure  bomb  assembly,  (b)  oven, 
and  (c)  automatic  control  and  temperature  recording  equipment.  The  sample 
was  contained  within  the  bomb,  which  was  placed  within  the  air  oven. 

The  temperature  of  the  sample  and  the  difference  between  it  and  the  oven 
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were  measured  with  a  thermocouple  system.  The  oven  temperature  was 
brought  slowly  to  a  value  at  which  self-heating  was  detected,  after  which 
automatic  control  minimized  the  temperature  difference.  Data  were  re¬ 
corded  until  the  experiment  terminated  when  the  sample  deflagrated  or 
exploded . 

Bombs  of  two  types  were  used;  schematic  drawings  of  these  are  shown 
in  Fig.  6a  (Type  i)  and  6b  (Type  II).  The  construction  material  was 
nickel,  chosen  because  of  the  possible  corrosive  properties  of  the  di- 
fluoroamino  propanes.  The  Type  I  bomb  utilized  the  three-thermocouple 
arrangement  of  Fig.  7a,  and  the  two  couples  in  the  bomb  were  retained  by 
modified  l/l6-in.  i.d.  Swagelok  male  connectors  (not  shown)  welded  in 
place.  The  lower  and  upper  bomb  casings  were  welded  after  the  bomb 
plug  was  in  place.  The  improved  Type  II  utilized  the  two-thermocouple 
arrangement  of  Fig.  7b,  and  the  couple  in  the  bomb  was  welded  to  the  lid. 
Type  II  is  an  improvement  because  of  its  greater  simplicity,  decreased 
mass,  and  more  nearly  spherically  symmetric  geometry. 

The  thermocouples  were  manufactured  by  Continental  Sensing,  Inc., 
Melrose  Park,  Delaware.  The  first  used  were  chromel-alumel  in  stainless 
steel  jackets.  Although  they  were  satisfactory,  copper-constantan  in 
Inconel  jackets  were  later  used  for  improved  sensitivity.  The  leads  for 
temperature  measurement  were  brought  to  a  Leeds  and  Northrup  Speedomax  H 
Azar  (adjustable  zero,  adjustable  range)  recorder  used  in  conjunction  with 
an  ice  bath.  Temperatures  were  measured  to  better  than  0.1°K.  The  leads 
for  differential  measurement  and  control  were  connected  to  a  Leeds  and 
Northrup  Model  9835  stabilized  d-c  microvolt  amplifier.  The  amplified 
voltage  was  recorded  and  used  for  proportional  control  of  the  nichrome 
heater.  For  this  purpose  a  Leeds  and  Northrup  Series  60  C.A.T.  system 
was  used  to  control  the  output  of  a  magnetic  amplifier.  Except  during 
periods  of  very  rapid  temperature  change,  temperature  differences  were 
detected  to  0.01°K  and  controlled  to  better  than  0.1°K. 

The  pressure  bomb  assembly  was  hung,  by  the  plug  extension  of  the 
thermocouple,  from  a  Transite  plate,  equipped  on  its  upper  surface  with 
a  spirally  wound  nichrome  wire  heater  and  a  centrifugal  blower.  Holes 
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WELDED 


MATERIAL:  COLD  DRAWN  NICKEL 
FINISH:  MACHINE 

(b)  TYPE  E 


RS  *  405*  *  IIP 


FIG  6  SCHEMATIC  DRAWINGS  OF  A  (Type  I)  AND 
B  (Type  II)  PRESSURE  BOMB  ASSEMBLIES 
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COPPER 

CONSTANTAN 


A-405I-I3R 


FIG.  7  THERMOCOUPLE  CONNECTIONS  USED  IN  PRESSURE  BOMB  ASSEMBLIES 

were  drilled  through  the  plate  at  an  angle  so  that  air  passing  through 
the  blower  could  pass  over  the  heater  to  the  bomb  without  the  necessity 
of  a  path  for  direct  radiation  between  heater  and  bomb.  In  addition  to 
the  couple  or  couples  within  the  bomb,  an  additional  couple  was  placed 
just  outside  of,  but  not  in  contact  with  the  bomb  (see  Fig.  7).  The 
assembly  was  set  into  a  Dewar  flask  for  which  the  Transite  plate  formed 
a  snug-fitting  lid,  with  the  bomb  at  the  center  of  the  flask,  which  was 
9  inches  in  diameter  and  10  inches  long.  The  flask  was  packed  in  vermiculite 
contained  in  a  heavy-walled  steel  pipe.  Addition  of  the  liquid  and  all 
subsequent  operations  were  performed  remotely  for  safety. 

The  experiment  was  carried  out  as  follows:  The  bomb  was  cleaned 
thoroughly  by  washing  successively  in  CH3C13,  HN03 ,  and  Ha0.  The  liquid 
sample  was  added  to  about  80$  of  the  bomb  capacity,  in  order  to  allow 
for  thermal  expansion  of  the  liquid;  this  amounted  to  2.1  ml  of  liquid 
at  room  temperature.  After  assembly,  heating  under  manual  control  over 
a  period  of  time,  tgjj,  of  about  5  hours  brought  the  sample  very  slowly 
to  a  temperature,  Tj,  at  which  there  was  barely  detectable  self-heating. 

At  that  point  automatic  control  of  the  temperature  difference  and  recording 
of  the  temperature  was  instituted  and  continued  until  the  experiment 
terminated. 
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2 .  Results 

Six  runs  were  made  with  1,2  DP,  and  three  with  2,2  DP.  Figure  8 
shows  typical  time-temperature  behavior,  the  data  being  that  of  run  ASH 
3-12;  the  figure  also  shows  the  plot  of  l/T  versus  In  dT/dt.  Figure  9 
shows  l/T  versus  In  dT/dt  for  four  runs  with  1,2  DP;  the  data  for  run 
ASH  13-12  resembled  the  other  runs  and  is  not  reproduced.  Figure  10 
shows  data  for  the  three  runs  with  2,2  DP.  Results  are  summarized  in 
Table  III. 

Table  III  shows  that  the  results  with  1,2  DP  are  not  reproducible  and 
yield  extremely  high  values  for  the  activation  energy.  An  apparent  two- 
step  behavior  occurred  in  four  of  the  six  runs.  It  is  shown  in  the  next 
section  that  results  of  this  type  could  result  from  an  autocatalytic 
reaction  in  which  a  reaction  with  low  activation  energy  occurs  after  de¬ 
velopment  of  sufficient  catalyst  concentration.  However,  the  occurrence 
of  surface  reactions  or  complicated  simultaneous  reactions  would  invali¬ 
date  Eqs .  (3-1)  and  (3-2)  and  consequently  the  analysis.  Because  of  the 
unreasonably  high  values  of  activation  energy  we  conclude  that  the  thermal 
decomposition  is  not  simple,  and  that  this  method  is  not  useful  for  de¬ 
termining  Eft .  Unless  a  new  approach  or  analysis  becomes  apparent,  no 
further  experiments  of  this  type  will  be  made  with  1,2  DP. 

In  contrast,  2,2  DP  behaved  in  accordance  with  expectations  based  on 
the  assumption  of  the  validity  of  the  method,  giving  an  average  value 
Ea  =  51.0  ±  2.5  kcal/mole.  Two  of  the  runs  were  made  In  Type  II  spherical 
bombs.  If  the  approximation  of  Eq.  (3-9)  Is  valid,  and  qv  can  be  calculal 
from  the  intercept,  Bq.  (3-8),  giving  log  qv  =  22.1  +  0.7. 

D .  Theoretical* 

Suppose  two  reactions  proceed  consecutively  without  loss  of  heat  to 
the  surroundings: 

X  L  +  Products  (3-11-1) 

X  +  L  Products  (3-11-2) 

*This  section  contributed  by  C.  M.  Ablow,  Manager  of  Mathematical  Sciences 
Department. 
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(o)  TYPICAL  TIME -TEMPERATURE  RECORD  IN 
ADIABATIC  SELF-HEATING  EXPERIMENT 
(RECORD  NUMBER  ASH  3-12). 


2150  2170  21.90  22.10  22.30 

T''x  I04 

(b)  PLOT  OF  In  dT/dt  v*  l/T 


•C  .  *OV  -M 


FIG.  8  (o)  TYPICAL  TIME-TEMPERATURE  RECORD 
IN  ADIABATIC  SELF-HEATING 
EXPERIMENT  (Record  No.  ASH  3-12) 

(b)  PLOT  OF  In  dT  dt  vs.  l/T 
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FIG.  10  PLOTS  OF  In  dT  dt  vs.  1  T  FOR  ADIABATIC  SELF-HEATING  EXPERIMENTS 
WITH  2,2  DP 
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It  is  shown  in  the  following  paragraphs  that  graphs  of  In  (dT/dt)  vs. 

( l/T)  have  the  qualitative  forms  sketched  in  Fig.  11,  where  Ej  and  E2  are 
the  activation  energies  of  the  two  reactions. 


FIG.  11  PLOTS  OF  In  dT/dt  vs.  l/T  FOR  CONSECUTIVE 
REACTIONS 


Equations  for  the  heat  and  mass  balance  are 

c^  =  qiVie  El/RT  +  q2v2  [L]e“E*/RT  (heat)  (3-12) 

|i=V,e-E'/»T  -  v2[L]e-E’/RT  (mass)  (3-13) 

where  [L]  is  the  concentration  of  reactant  L,  and  subscripts  1  and  2  in¬ 
dicate  the  reactions  of  Eqs.  (3-11-1)  and  (3-11-2)  respectively. 


From  Eqs.  (3-12)  and  (3-13)  one  may  derive  the  differential  equation 
valid  in  the  (T,[l])  phase  plane: 


d  [l] 

dT 


(3-14) 


or 


d[L] 

dT 


-AE/RT  ,  , 

ye  '  -  [L] 

p  -AE/RT  TT 
qve  '  +  ILJ 
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where 

AE  —  E 1  —  E2 , 

q  =  qi/q2. 

V  =  Vi/v2 . 

and  p  =  c/q2  . 


Equation  (3-15)  gives  the  slope  of  the  graph  of  [h]  vs.  T  at  any 
point  in  the  ( [L] ,T)-plane .  In  particular  the  curves  connecting  points 
at  which  that  slope  is  constant  may  be  plotted.  If  N  is  the  numerator 
of  the  right  hand  side  of  Eq .  (3-15) 

N  =  p  [Ve_AE^RT  -  [L]}  (3-16) 


then  the  curve  N  =  0  connects  points  with  zero  slope, 
is  the  denominator 


D  = 


-ae/rt 

qve  + 


[L] 


Similarly  if  D 

(3-17) 


then  the  curve  D  =  0  connects  points  with  infinite  slope.  The  whole  set 
of  solutions  of  [ L ]  vs.  T  for  various  initial  conditions  may  now  be 
plotted.  The  times  corresponding  to  points  on  a  particular  ([l],T)  solu¬ 
tion  curve  are  obtained  by  integrating  the  mass  balance  Eq .  (3-13), 

Three  cases  are  to  be  distinguished,  AE  <  0,  AE  =  0,  and  AE  >  0 . 
These  cases  are  shown  graphically  in  Figs.  12,  13  and  14. 


Case  I:  AE  <  0 .  The  solution  curve  is  similar  to  BCFG  in  Fig.  12. 
At  time  zero,  point  B,  the  temperature  is  low  and  the  concentration  of  L 
is  zero.  As  time  goes  on  the  temperature  and  concentration  [L]  both  rise 
along  BC .  The  concentration  rises  to  a  maximum  at  F  and  then  falls  a 
little  toward  G  as  the  temperature  continues  to  rise. 


The  time  rate  of  temperature  rise  can  be  found  from  the  heat  balance 


Eq.  (3-12): 


or 


dT _ £  -E2/RT 

dt  -  v* 


(qve 


-ae/rt 


+  [L]} 


(3-18) 


=  _E  e-E*/RTD 


dT 


dt  v2 


(3-19) 
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FIG.  12  PLOTS  OF  CONCENTRA7  ION  vs. 
TEMPERATURE  FOR 
CONSECUTIVE  REACTIONS  (AE  <  0) 


FIG.  13  PLOTS  OF  CONCENTRATION  vs. 
TEMPERATURE  FOR 
CONSECUTIVE  REACTIONS  (AE  =  0) 


FIG.  14  PLOTS  OF  CONCENTRATION  vs. 
TEMPERATURE  FOR 
CONSECUTIVE  REACTIONS  (AE  >  0) 
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Near  time  zero,  [l]  is  near  zero  and  the  temperature  is  low.  Both  these 
properties  tend  to  make  the  exponential  term  in  the  braces,  Eq .  (3-18), 
dominate  the  concentration  [L].  In  approximation,  then,  near  time  zero 

dT  =  gqv  e“Ei/RT.  (3-20) 

dt  V2 

The  bracket  term  is  D  of  Eq .  (3-17).  Curves  of  constant  D  are  parallel 
in  the  sketch  to  the  curve  D  =  0  shown.  Thus  at  later  times,  along  F  G  in 
the  sketch,  D  is  approximately  constant.  At  late  times  Eq .  (3-19),  with 
D  constant  is  a  good  approximation. 

In  summary,  a  plot  of  ln(dT/dt)  vs.  (l/T)  would  have  slope  near 
(-E i/R)  at  early  times  and  near  (-Ea/R)  for  late  times  in  this  case  where 
E2  >  . 

Case  II:  AE  =  0.  Both  [L]  and  T  continually  rise  during  the  reaction. 
The  temperature  rises  without  limit  while  the  concentration  of  reactant 
L  rises  toward  an  upper  bound.  Equation  (3-18)  shows  that  near  time  zero, 
while  [ L ]  is  much  smaller  than  qv ,  Eq.  (3-19)  with  D  constant  is  a  fair 
approximation.  At  late  times  Fig.  13  shows  that  [l]  becomes  nearly 
constant  so  that  again  Eq .  (3-19)  is  approximately  correct  with  D  a 
larger  constant. 

In  this  case,  a  plot  of  ln(dT/dt)  vs.  (l/T)  has  slope  (-Ej/R)  both 
for  early  and  late  times,  with  an  intermediate  region  of  transition  at 
greater  slope  from  the  lower  of  two  parallel  lines  to  the  upper. 

Case  III:  AE  >  0.  Near  time  zero,  [l]  in  the  bracket  in  Eq .  (3-18) 

is  negligible  and  Eq,  (3-20 )  is  valid.  At  late  times  Fig.  14  shows  that 

the  two  terms  of  the  bracket  approach  constants  so  that  Eq .  (3-19)  with 

D  constant  is  an  approximate  solution.  In  this  case  the  slope  of  a  plot 

of  ln(dT/dt)  vs.  (l/T)  drops  with  time  from  (-Ej/R)  to  (-E2/R). 
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E.  Future  Plans 


Additional  experiments  will  be  carried  out  with  2,2  DP  in  pressure 
bombs  with  thinner  walls  in  order  to  demonstrate  the  validity  of  Eq. 
(3-10)  and  to  develop  a  better  model  for  calculating  reliable  values  of 
Ea  and  qv . 

Should  these  experiments  confirm  the  value  of  the  self -heating 
method,  studies  w:ll  be  initiated  on  additional  compounds. 
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XV  KINETICS  OF  THERMAL  DECOMPOSITION 
(Theodore  Mill) 


A .  Summary 


The  objective  of  this  phase  is  to  acquire  a  knowledge  of  the  nature 
and  rate  of  the  rate-controlling  step  in  decompositions,  and  of  the  nature 
of  intermediates,  of  transition  states,  and  of  the  ultimate  products  in 
order  to  devise  a  description  of  the  mechanism  of  the  reaction. 

Kinetic  and  mechanism  experiments  were  carried  out  with  1,2  DP  as 
a  neat  liquid,  as  a  dilute  solution  in  nitrobenzene,  and  in  the  vapor 
phase.  A  few  experiments  were  carried  out  on  the  decomposition  of  2,2  DP 
in  the  vapor  phase.  The  two  compounds  apparently  differ  but  little  in 
their  kinetic  behavior,  but  have  significantly  different  rates. 

The  decomposition  of  1,2  DP  in  nitrobenzene  solution  in  the  tempera¬ 
ture  range  176-187°C  is  described  by  the  equation 

rate  =  dx/dt  =  kx (A  -  x)  +  k2 (A  -  x)(x)2 

where  (A  -  x)  is  the  instantaneous  concentration  of  1,2  DP.  Thus  the 
reaction  is  autocatalytic  over  60-70$  of  the  decomposition.  Kinetic 
parameters  for  the  solution  decomposition  are 

kj  =  3.7  x  1014  exp[  (-40.2  t  2)/RT  sec*"1 

k2  =  1.8  x  1014  exp [(-35. 7  ±  2)/RT]  liter2mol“2 sec'1 

where  the  units  of  Ea  are  kcal/mole.  The  vapor  phase  decomposition  dis¬ 
plays  marked  surface  effects  in  glass  and  steel  and  is  also  autocatalytic. 

2,2  DP  decomposes  at  about  one-tenth  the  rate  of  1,2  DP.  The  de¬ 
composition  displays  autocatalysis  both  in  solution  and  in  the  vapor  phase; 
there  are  marked  surface  effects  in  the  vapor  phase. 
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B ,  Mechanism  and  Kinetics  of  the  Thermal  Decomposition  of  1,2  DP 


1 .  Kinetics  Studies 

a .  Decomposition  of  Pure  Liquid 

1,2  DP  containing  1  to  2  percent  methylene  chloride,  was  stored 
on  a  vacuum  line  in  small  tubes  inside  polyurethane  "totes."  Twenty  to 
thirty  mg  samples  were  frozen  into  tared  capillary  tubes  on  the  vacuum 
line.  The  tubes  were  sealed  with  a  flame  and  immersed  for  varying  in¬ 
tervals  of  time  in  a  bath  at  153°C.  Contents  of  each  capillary  were 
analyzed  by  gas  chromatography.  Several  volatile  products  were  found, 
but  none  was  present  in  an  amount  equal  to  the  loss  of  1,2  DP.  On  long 
heating,  the  only  product  was  a  dark  nonvolatile  residue.  The  results  of 
these  experiments  are  given  in  Table  IV.  The  first  half-life  is  about  20 
hours  and  the  rate  markedly  accelerates  thereafter.  Since  the  polymeric 
residue  separates  from  the  starting  material  at  an  early  stage  (12-20  hr), 
the  actual  concentration  of  1,2  DP  may  not  change  appreciably,  although 
the  total  amount  present  does  change. 

Table  IV 

DECOMPOSITION  OF  PURE  LIQUID  1,2  DP  AT  153°C 


Tube 

No. 

Time 

(min) 

Initial 

Weight 

(mg) 

Recovered 

Reactant 

(">g) 

Percent 

Decomposition 

Residue 

(mg) 

n 

1260 

34.7 

22.9 

34 

_ 

1440 

33.7 

12.7 

62 

4.6 

1860 

30.0 

2.5 

94 

7.3 

m I 

2535 

21.6 

0.05 

99+ 

22.1 

There  are  a  number  of  experimental  difficulties  associated  with 
a  kinetic  study  of  the  pure  liquid.  The  low  boiling  point  (81°C)  and 
the  low  rate  of  decomposition  require  the  use  of  high  temperature,  and 
of  pressures  in  excess  of  seven  atmospheres,  to  provide  conveniently  fast 
rates  of  decomposition;  thus  container  and  sampling  problems  are  difficult. 
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Furthermore,  the  uncertainty  associated  with  the  concentration  of  the 
neat  liquid  during  decomposition  and  the  invariance  of  the  initial  con¬ 
centration  place  limitations  on  the  possible  kinds  of  kinetic  experiments. 
We  therefore  decided  to  examine  the  decomposition  in  dilute  solution, 
in  the  expectation  of  returning  later  to  the  problem  of  neat  liquid  de¬ 
composition  by  direct  experimentation  and  by  extrapolation  of  the  solution 
study  results , 

b .  Decomposition  of  Solution  in  Nitrobenzene 

Nitrobenzene  was  chosen  as  the  solvent  for  the  solutions  studies 
because  of  its  high  boiling  point  (211°C),  inertness,  good  solvent  pro¬ 
perties,  and  ease  of  purification.  1,2  DP  in  nitrobenzene  in  concentra¬ 
tions  between  0.7  and  0.1  molar  was  found  to  decompose  at  a  convenient 
rate  for  study  in  the  temperature  range  175-190°C.  Samples  were  vacuum 
distilled  into  tared  vessels  containing  known  amounts  of  vacuum-distilled 
nitrobenzene.  The  gross  features  of  the  reaction  were  the  same  as  those 
observed  in  the  decomposition  of  the  neat  liquid.  A  polymeric  phase 
separated  during  the  decomposition;  the  same  intermediates  formed  and  dis¬ 
appeared,  and  the  rate  of  decomposition  accelerated  after  about  one 
quarter-life.  A  gas  chromatogram  of  a  partially  decomposed  solution  is 
shown  in  Fig,  15. 

Rate  measurements  were  made  in  two  ways.  For  both  methods  a 
standard  solution  of  1,2  DP  in  nitrobenzene  was  used  to  standardize  the 
samples.  Reproducibility  of  the  standard  was  about  2  percent  from  run  to 
run  during  one  day.  In  the  first  method,  10  pi iter  samples  of  solution 
were  withdrawn  from  1.5  ml  of  bulk  solution  in  glass  through  a  serum  cap 
and  Teflon-glass  stopcock.  Aliquots  were  injected  directly  into  the  gas 
chromatograph  and  the  amounts  of  1,2  DP  were  estimated  by  comparison  with 
the  chromatogram  of  a  standard  sample.  These  data  were  often  not  repro¬ 
ducible,  because  of  evaporative  loss  of  the  1,2  DP  through  the  stopcock 
and  septum.  This  loss  was  demonstrated  by  following  the  change  in  con¬ 
centration  of  pure  n-octane  dissolved  in  nitrobenzene.  Results  of  several 
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FIG.  15  CHROMATOGRAM  OF  DECOMPOSITION 
PRODUCTS  OF  NITROBENZENE 
SOLUTION  OF  1,2  DP  OBTAINED  ON  A 
5-FT  DINONYLPHTHALATE  COLUMN 
AT  70  C  AND  30  ml/min  FLOW  RATE 


runs,  made  at  176°C  under  conditions  for  which  no  loss  of  octane  was 
noted,  are  given  in  Fig.  16. 

The  second  method  involved  the  preparation  of  a  number  of 
sealed  glass  tubes,  each  containing  50  joliter  of  solution.  Pairs  of  tubes 
were  removed  from  the  temperature  bath  at  selected  time  intervals  and 
the  contents  analyzed  by  gas  chromatograph  .  Initial  values  for  dupli¬ 
cate  samples  agreed  closely  but  divergence  in  results  increased  markedly 
as  the  reaction  rate  accelerated.  Results  of  runs  made  at  187°C  are 
given  in  Fig.  17. 
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FIG.  16  MOLAR  CONCENTRATION  OF  1,2  DP  vs.  TIME 
FOR  NITROBENZENE  SOLUTION  OF  1,2  DP 
AT  176  C  (Aliquot  Method) 


TIME - hr* 


FIG.  17  MOLAR  CONCENTRATION  OF  1,2  DP  vs.  TIME 
FOR  NITROBENZENE  SOLUTION  OF  1,2  DP 
AT  187  C  (Tored  Somple  Method) 
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c  *  Decomposition  of  Vapor 

Experiments  were  carried  out  in  calibrated  50-ml  glass  and 
steel  bulbs.  Known  amounts  of  1,2  DP  were  added  to  each  bulb  by  freezing 
into  the  bulb  quantities  of  material  from  a  measured  pressure  of  the  vapor 
in  a  vacuum  manifold  of  known  volume.  Bulbs  were  removed  from  the  oil 
bath  at  intervals,  the  contents  frozen  in  liquid  nitrogen,  and  noncon¬ 
densables  estimated  by  pressure  rise  on  the  vacuum  line.  Samples  of  the 
noncondensable  products  were  analyzed  by  mass  spectrometry;  only  trace 
amounts  of  H2  and  N2  were  found.  Condensables  were  transferred  on  the 
line  into  nitrobenzene  and  analyzed  by  gas  chromatography. 

The  volatile  products  were  identical  to  those  observed  in  liquid 
phase  and  solution  decompositions;  there  was  a  significant  increase  in 
the  amount  of  peak  2  of  Fig.  15.  A  brown  polymer  was  deposited  on  the 
walls  of  the  bulbs.  The  vapor  phase  decomposition  rate  was  about  one- 
tenth  that  of  the  nitrobenzene  solution.  There  was  a  marked  surface 
effect,  which  prevented  obtention  of  reproducible  results  either  in  glass 
or  steel.  The  surface  effect  takes  the  form  of  an  increase  in  decomposi¬ 
tion  rate  in  used  steel  and  a  retardation  in  used  glass.  Probably  the 
reaction  takes  place  on  the  walls.  Table  V  gives  approximate  half-lives, 


Table  V 

VAPOR  PHASE  DECOMPOSITION  OF  1,2  DP 


T 

(°c) 

Concentration  x  103 
(mole/l ) 

Container 

tl/2 

(hrs) 

0.64 

New  glass 

-  80 

0.64 

Used  glassa 

»  80 

1 

0.64 

Used  glass** 

-  55 

mSm 

1.1 

Used  stainless  steel 

<  20 

aBulbs  were  used  for  one  previous  run  at  177°C 
^Bulbs  were  used  for  two  previous  runs  at  176°C. 
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2.  Mechanism  Studies 


a .  Catalysis  by  Water 

Added  water  was  found  to  cause  a  marked  increase  in  decomposi¬ 
tion  rate  of  the  pure  liquid  at  153°C  but  no  quantitative  data  were  ob¬ 
tained.  The  product  was  a  dark  polymer. 


The  effect  of  water  on  the  rate  of  decomposition  of  1,2  DP  in 
nitrobenzene  solutions  was  studied  at  176°C.  One  solution  studied  was 
0.227  M  with  respect  to  water  and  0.655  M  in  1,2  DP;  the  other  was  0.019 
M  in  water  and  0.576  M  in  1,2  DP.  This  latter  small  and  uncertain  amount 
of  water  was  present  originally  in  the  1,2  DP  as  obtained  from  the  dis¬ 
tillation.  Rate  measurements  were  made  using  sealed  tubes  with  duplicate 
samples  for  each  point,  with  the  results  shown  in  Fig.  18.  There  is  a 
marked  increase  in  rate  for  the  "wet”  (0.227  M)  solution  compared  to  that 
for  the  ”dryM  (0.019  M)  solution.  The  apparent  first-order  rate  constant 
for  the  wet  solution  decomposition  was  estimated  to  be  4  x  10“5  sec1"1  , 
which  is  more  than  three  times  the  rate  for  the  runs  reported  in  Fig.  16. 
The  apparent  first  order  rate  constant  for  the  dry  run  is  1.6  x  10~5  sec"1 
or  about  twice  the  value  of  Fig.  16,  an  increase  of  rate  which  may  be  due 
to  the  presence  of  water. 

b .  Effect  of  Inhibitors 

Two  kinds  of  experiments  were  carried  out  to  test  for  the 
presence  of  a  kinetic  chain  (free-radical  or  ionic)  in  the  decomposition. 

A  0.67  M  solution  of  1,2  DP  in  nitrobenzene  containing  5  mole  percent  of 
2 ,6-di-t-butyl-4-methyl phenol ,  a  well-known  free  radical  scavenger,16 
was  heated  for  5  hours  at  163°C,  The  solution  became  very  dark  and  66 
percent  of  1,2  DP  disappeared.  AIM  solution  without  inhibitor  decomposed 
more  slowly,  with  70  percent  decomposition  in  10  hours.  A  nitrobenzene 
solution  of  the  inhibitor  alone  gave  no  color  change  on  heating  for  twenty 
hours . 

16G.  S.  Hammond,  J.  N.  Ser,  and  C,  E.  Boozer,  J.  Amer.  Chem.  Soc.,  77, 

3244  (1955) 
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FIG.  18  FIRST  ORDER  PLOTS  FOR  THERMAL  DECOMPOSITION 

OF  1,2  DP  IN  DRY  (0.019  M  water)  AND  WET  (0.227  M  water) 
NITROBENZENE  SOLUTIONS 
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A  0.5  M  solution  of  1,2  DP  in  a-methyl styrene ,  an  efficient 
radical  and  proton  trap,17  was  sealed  in  an  NMR  sample  tube.  After  heat¬ 
ing  at  163°C  for  5  hours,  the  solution  was  very  dark.  An  F19  NMR  spectrum 
of  the  contents  showed  an  absence  of  bands  due  to  1,2  DP18,  but  the  presence 
of  two  weak  bands  well  upfield  (128  and  148  ppm  from  CC13F)  from  the  NF2 
bands  of  1,2  DP.  The  position  of  these  new  bands  suggests  they  are  due 
to  carbon-bonded  fluorine.  Other  F19  nuclei  may  have  escaped  detection 
in  the  spectrometer  because  of  band  broadening  in  the  tarry  product  or 
because  of  loss  of  volatile  SiF4  into  the  void  above  the  liquid. 

The  acceleration  in  rate  found  with  both  the  di-butyl  methyl- 
phenol  and  the  a-methylstyrene  seems  to  rule  out  the  possibility  of  a 
free-radical  chain  decomposition. 

c ‘  Container  Effects 

Although  wall  effects  in  solution  kinetics  are  unusual,  this 
system  may  involve  formation  of  hydrogen  fluoride,  and  its  interaction 
with  glass  or  metal  oxides  could  introduce  spurious  kinetic  effects.  To 
assess  the  importance  of  wall  effects,  0.3-0. 5  M  solutions  of  1,2  DP 
in  nitrobenzene  were  heated  in  stainless  steel  and  glass  vessels  and  in 
a  glass  tube  packed  with  glass  wool .  Rates  were  followed  by  removing 
aliquots  at  various  intervals.  No  great  difference  in  initial  rates 
was  found;  the  packed  tube  rate  was  twice  that  of  the  other  glass  or  steel 
vessels.  Differences  between  steel  and  glass  were  no  greater  than  among 
different  runs  in  glass. 

Several  runs  were  made  using  sealed  steel  and  brass  vessels. 

The  product  of  decomposition  was  grossly  the  same  in  all  runs,  a  dark 
polymer  with  the  odor  of  hydrogen  fluoride. 


17F.  R.  Mayo  and  C.  Walling,  Chem  Revs.  46,  191  (1950) 

18F19  NMR  spectra  of  difluoroamino  compounds  are  reported  In  Rohm  and 
Haas  Quarterly  Progress  Report  P-60-10,  May  31,  1960,  p.  24  (Conf . ) ; 
all  compounds  exhibit  bands  at  -4000  to  -5000  cps  from  CF3C00H. 

Z  has  bands  at  -2134  cps  and  -3219  cps  from  CC13F  (-37.8  and  -57.1 
ppm  from  CC13F). 
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d .  Decomposition  Products 

Several  intermediate,  volatile  products  were  observed  in  gas 
chromatograms.  At  no  time  did  the  total  mass  of  these  products  appear  to 
account  for  the  loss  of  1,2  DP  and  on  prolonged  heating  they  disappeared, 
an  indication  of  their  intermediate  character.  (Water  cannot  be  detected 
with  the  chromatographic  column  used  for  these  analyses.)  Two  peaks  in 
the  chromatogram  were  present  in  significant  amounts.  These  are  shown  in 
Fig.  15  labeled  3  and  4.  Methylene  chloride,  the  solvent  from  which 
1,2  DP  is  isolated,  was  normally  present  to  the  extent  of  1  percent  in 
distilled  1,2  DP.  The  compound,  free  of  CH2C12,  behaves  in  an  identical 
way  to  the  other  samples.  Peaks  3  and  4  were  isolated  by  gas  chroma- 
tography  and  their  mass  spectra  obtained.  Peak  3  has  the  cracking  pattern 
of  acetonitrile  (CH3CN)  and  a  retention  time  identical  to  CH3CN  in  the 
chromatograph.  Peak  4  displays  a  significant  peak  at  m/e  =  86  in  the 
mass  spectrum.  We  were  unable  to  clear  the  spectrum  of  possible  impurities 
(l,2  DP,  water,  acetone,  air)  to  give  residual  pattern  in  reasonable  ac¬ 
cord  with  structures  having  a  parent  mass  of  86.  One  such  reasonable 
structure  is  CH3^CN. 

Several  samples  of  polymer  have  been  isolated  from  the  decom¬ 
positions  of  1,2  DP  in  liquid  solution,  and  vapor  phase.  None  of  the 
analyses  is  satisfactory  for  CHN  balance;  no  fluorine  analysis  has  been 
obtained.  Water  very  probably  contributes  some  oxygen  to  the  polymer. 

Most  samples  of  polymer  are  insoluble  in  acetone,  suggestive  of  a  highly 
condensed  system  with  few  functional  groups  such  as  carbonyl  or  hydroxyl. 

e .  Solvent  Effects 

In  one  experiment  hexadecane  was  substituted  for  nitrobenzene. 

A  sealed  NMR  tube  was  used  as  the  container  and  the  F19  spectrum  of  a 
sample  heated  18  hours  at  163°C  showed  little  change  in  the  intensity  of 
the  two  fluorine  bands,  compared  to  the  spectrum  of  the  unheated  solution. 
Both  spectra  were  obtained  at  about  the  same  time.  A  similar  solution 
in  nitrobenzene  had  half  decomposed  in  6  to  7  hours. 
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f  .  Reaction  with  Alkoxide  Ion 

One  mmole  of  1,2  DP  dissolved  in  n-butanol  (chosen  for  conven¬ 
ience  in  the  gas  chromatographic  analysis)  was  treated  with  4  mmoles  of 
sodium  butoxide  in  butanol  over  a  12-hour  period  at  25°C.  Samples  were 
chromatographed  at  intervals.  1,2  DP  was  very  rapidly  consumed  after 
each  addition  of  base,  and  product  peaks  which  appeared  in  the  chromato¬ 
grams  corresponded  to  only  a  small  fraction  of  the  consumed  starting 
material .  These  products  are  identical  with  those  observed  in  thermal 
decomposition  experiments  (Fig.  15)  with  respect  to  both  retention  times 
and  peak  shapes. 

3 .  Discussion 

a  .  Rate  Equation  and  Rate  Constants  for  Thermal 
Decomposition  in  Nitrobenzene  Solution 

The  results  shown  in  Figs.  12  and  13  indicate  that  the  decom¬ 
position  of  1,2  DP  in  nitrobenzene  is  initially  first  order,  with  ac¬ 
cumulation  of  an  intermediate  or  intermediates  producing  autocatalysis. 
Figure  16  gives  composite  results  of  several  runs  at  176°C,  as  plots  of 
log  concentration  vs.  time.  The  data  yield  good  straight  lines  for  20-30 
percent  of  reaction.  Over  an  almost  five-fold  change  in  initial  concen¬ 
tration  from  0.56  to  0.14  moles/liter  the  initial  slopes  are  identical 
within  experimental  error  and  give  a  value  kx  =  1.0  x  10“5  sec-1  at  176°C. 
At  187°C  the  plot  of  log  concentration  versus  time  gave  similar  results 
with  greater  scatter  (Fig.  17).  By  a  least  squares  calculation  for 
straight  line  fits  of  the  points  over  the  first  20-30  percent  of  decom¬ 
position,  the  value  of  kj  was  found  to  be  3.0  x  10“5  sec”1  at  187°C. 

Assign  A  as  the  original  concentration  of  1,2  DP,  and  x  the 
amount  which  has  been  converted  into  product  at  any  given  time.  Then 
if  one  assumes  that  acceleration  is  caused  either  by  a  primary  product 
or  by  a  secondary  product  produced  rapidly  compared  to  the  rate-controlling 
step,  the  rate  is  described  by  the  equation 

rate  =  dx/dt  =  kl (A  -  x)  +  k2(A  -  x)(x)y.  (4-1) 
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nkj , 


Eq .  (4-l)  becomes 


(4-2) 

(4-3) 


dx/dt  -  (A  -  x)ki(l  +  nxy) 

If  we  define  an  average  rate  E  =  Ax/At,  then  Eq .  (4-3)  becomes 

[R/kj (A  -  x ) ]  -  1  =  nxy  (4-4) 


Table  VI  gives  calculated  values  of  x,  R,  and  kt (A  -  x)  for 
run  at  176°C  with  a  0.461  M  solution  in  nitrobenzene,  assuming  kx  =  1 
10*"5  sec-1.  A  plot  of  these  data  in  terms  of  log  R/kj  (A  -  x)  vs.  log 
shown  in  Fig.  19.  After  7  hours  the  points  fall  on  a  straight  line. 


Table  VI 

RESULTS  OF  KINETICS  EXPERIMENT  WITH  0.461  M  1,2  DP 
IN  NITROBENZENE  AT  176°C 


Time 
(hrs ) 

x  •  10~2 
(moles/liter) 

R  •  102 

(moles/liter/hr) 

kt(A  -  x)  *  102 
(moles/liter/hr ) 

0 

0 

1.4 

1.6 

1 

1.4 

1  .5 

1.56 

2 

2.9 

1.6 

1.51 

3 

4  .5 

1.5 

1.46 

4 

6.0 

1.6 

1  .40 

5 

7.6 

1 .8 

1.35 

6 

9.4 

2.4 

1 .28 

7 

10.0 

3.1 

1.20 

8 

14.9 

3.8 

1.09 

9 

18.7 

4.3 

0.959 

10 

23.0 

4.5 

0.809 

11 

27.5 

4.4 

0.651 

12 

31.9 

4.3 

0.497 

13 

36.2 

0.347 
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FIG.  19  LOG-LOG  PLOT  OF  lR^A-x);  -1  *  nxy  FOR  0.416  M  SOLUTION 
OF  1,2  DP  IN  NITROBENZENE  AT  176  C 


R  and  x  change  most  slowly  during  the  first  7  hours,  the  changes  in  x 
are  subject  to  the  greatest  experimental  error  in  that  period.  For 
this  reason,  the  early  points  were  ignored  in  plotting  the  straight  line. 
The  slope  gave  a  value  of  y  =  2.07.  When  the  line  was  adjusted  to 
y  =  2.0,  the  intercept  gave  n  =  83 .  A  similar  analysis  was  carried  out 
for  the  runs  at  187°C  and  for  a  value  y  =  2,  n  =  64  was  obtained.  By 
Eq.  (4-2)  ka  is  found  to  be  8.2  x  10“4  liter2/mole“a/sec"1 .  Figure  20 
shows  a  comparison  of  the  curves  by  plotting  (a)  R  from  Table  VI  versus 
time,  shown  as  circles,  and  (b)  Eq.  (4-3)  for  n  s  83,  y  =  2,  shown  as 
triangles.  The  correlation  is  reasonably  good.  Thus  the  following 
equation  describes  the  over-all  decomposition  at  176°C  in  nitrobenzene: 
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FIG.  20  PLOTS  OF  RATE  v*.  TIME  FOR  0.46  M  1,2  DP 

IN  NITROBENZENE  AT  176  C  -  EXPERIMENTAL 
DATA  AS  SHOWN  IN  TABLE  VI  -  EQ.  (4.3)  FOR 
n  ■  80,  y  •  2 
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dx/dt  =  1.0  x  10*"5  (A  -  x)(l  +  83x2)  moles/liter'"1/sec“1  (4-5) 

Similar  treatment  of  the  run  at  187°C  gave 

dx/dt  =  3.0  x  10“5  (A  -  x)(l  +  64x2  )  moles/liter^/sec-1  .  (4-6) 

Arrhenius  parameters,  v  (frequency  factor)  and  Ea  (activation 
energy)  were  calculated  from  the  values  of  kj  and  k2  at  176  and  187°C. 

They  are  given  in  Table  VII. 

Table  VII 

KINETIC  PARAMETERS  FOR  THERMAL  DECOMPOSITION 
OF  1,2  DP  IN  NITROBENZENE 


Rate 

T 

Ea 

V 

Constant 

(°K) 

(kcal/  mole ) 

( sec”1 ) 

ki 

449 

460 

40.2  +  2 

3.7  x  1014 

k2 

449 

460 

35.7  *  2 

X .8  x  1014 

b .  Decomposition  Mechanism  in  Nitrobenzene  Solution 

In  principle  several  unimolecular  decomposition  pathways  are 
available  to  1,2  DP.  Thermochemical  considerations  suggest  that  either 
an  ionic  (polar)  dissociation  involving  loss  of  hydrogen  fluoride  or  a 
free  radical  cleavage  of  a  C-N  bond  are  the  most  probable  modes  of  de¬ 
composition.  The  bulk  of  the  experimental  evidence  supports  a  mechanism 
involving  a  rate-controlling  loss  of  hydrogen  fluoride. 

The  diminution  in  rate  in  going  from  the  vapor  phase  (where 
the  reaction  probably  occurs  on  the  walls)  to  hexadecane  and t o  nitro¬ 
benzene  solutions  is  indicative  of  the  formation  of  a  polar  transition 
state,  or  increased  solvation  of  products  compared  to  reactants,  or  both.19 
An  alternative  explanation  that  there  is  change  in  mechanism  accompanying 
a  change  in  phase  or  solvent  is  unlikely  in  view  of  the  similarity  of  the 
products  formed  and  the  over-all  kinetic  behavior  exhibited  under  these 
different  conditions.  The  observation  that  butoxide  ion  gives  the  same 

l9A.  A.  Frost  and  R.  G.  Pearson,  Kinetics  and  Mechanism,  John  Wiley  and 
Son,  Inc.,  New  York,  1953,  p.  117. 

46 


CONFIDENTIAL 


CONFIDENTIAL 


kinds  of  products  at  25°C  as  are  found  in  the  thermal  decomposition  pro¬ 
vides  strong  support  for  a  mechanism  involving  loss  of  hydrogen  fluoride. 

Two  of  the  intermediate  products  observed  in  all  four  kinds 
of  decompositions  (vapor  phase,  pure  liquid,  nitrobenzene  solution,  and 
butoxide  ion  butanol  solvent)  have  been  identified  as  acetonitrile 
(CH3CN )  and  (tentatively)  2-(fluoroimino)propionitrile  (CH3C(=NF)CN) . 

Both  products  result  from  loss  of  hydrogen  fluoride  from  1,2  DP: 


ch3c-ch2nf2 

nf2 


HF  +  CH3CCH2NF2 
NF 


CH3fCH2NF2 

NF 


CH3CN  +  FCH2NF2 
CH3C(NF)CN  +  2HF 


Neither  free-radical  inhibitor  was  found  to  inhibit  the  decom¬ 
position  of  1,2  DP  in  solution,  a  result  in  keeping  with  the  suggested 
mechanism.  Unexpected  acceleration  of  decomposition  found  with  both  in¬ 
hibitors  probably  is  due  to  the  intervention  of  a  fast  base-catalyzed 
elimination  of  HF  by  the  conjugate  anion  of  di-t-butyl-methylphenol ,  in 
one  case,  and  by  a  positive  solvation  effect  in  the  case  of  a-methyl- 
styrene.  a-Methyl styrene  could  also  have  some  nucleophilic  character 
towards  1,2  DP;  a  reasonable  transition  state  structure  can  be  written: 


The  over-all  rate  equation  for  the  decomposition  (Eq.  4-5) 
predicts  a  second-order  dependence  on  the  catalytic  species.  Since  a 
termolecular  collision  is  improbable,  the  equation  implies  that  the 
catalytic  species  may  be  a  secondary  product  derived  from  a  first  formed 
product,  possibly  HF,  present  in  amount  x  at  time  t: 

A  ->  A  -  x  +  xB 
2B  C 

A  +  C  -»  A-x  +  xB  +  C 

In  this  scheme,  the  concentration  of  C  is  a  function  of  Ba ,  that  is  x* . 
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Water  is  known  to  catalyze  this  decomposition,  a  fact  which  is 
consistent  with  the  proposed  mechanism.  Probably  this  catalysis  is 
associated  with  water1 s  nucleophilic  character  rather  than  with  any  al¬ 
teration  of  the  solvating  properties  of  the  solvent.  Since  water  is  a 
secondary  product  in  the  reaction,  produced  by  the  reaction 

2HF  +  Si02  _>  H20  +  F2SiO, 

it  is  of  interest  to  inquire  whether  water  could  be  the  product  responsible 
for  autocatalysis.  The  experiment  in  which  water  was  added  initially  to 
the  nitrobenzene  solution  provides  a  way  of  testing  this  hypothesis. 

The  concentration  of  water  is  substituted  into  Eq .  (4-1 )  to  the 
first  power 

rate  =  dx/dt  s  k*(A  -  x)  +  k2 (A  -  x)(H20).  (4-7) 

Solving  for  the  initial  rate,  i.e,,  when  (A  -  x)  =  (A),  in  a  nitrobenzene 
solution  containing  0.655  M  1,2  DP  and  0.277  M  H20,  one  obtains  a  value 
for  the  initial  rate  of  1.56  x  10“4  mole/l iter“1/sec"1 .  This  is  almost 
four  times  the  measured  value  of  4  x  10“5  mole-1 iter-1 /sec  actually  found 
for  the  water-catalyzed  reaction,  and  the  difference  is  well  beyond  the 
experimental  error.  One  therefore  concludes  that  water  is  not  a  suf¬ 
ficiently  strong  nucleophile  to  account  for  the  observed  results. 

The  marked  surface  effects  found  in  vapor-phase  decompositions 
must  reflect  the  fact  that  the  reaction  occurs  largely  at  the  walls.  In 
steel  bulbs,  HF  produces  ferric  fluoride,  a  Lewis  acid  which  presumably 
can  catalyze  the  initial  decomposition  by  associating  with  a  nitrogen- 
bound  fluorine 

5+  &" 

N--F-- FeFs 

Glass  bulbs  exhibit  an  inhibiting  effect  with  continued  use.  It  may  be 
that  deposition  of  an  organic  film  on  the  glass  simply  covers  up  active 
sites  already  present  on  the  walls.  Surface  effects  have  also  been 
found  in  the  vapor-phase  decomposition  of  IBA  by  Booman.20  He  recently 

2 °K .  Booman.  Private  communication  to  M.  E.  Hill,  March  15,  1963. 
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reported  new  evidence  indicating  that  IBA  decomposes  by  loss  of  HF  rather 
than  by  a  radical  cleavage,  and  that  the  rate  is  accelerated  in  steel 
vessels  (FeF3  catalysis)  but  inhibited  in  brass  vessels. 

c .  Speculation  on  Mechanisms  and  Rates  in  Pure  Liquid 
and  at  High  Temperature 


It  is  of  interest  to  consider  the  applicability  of  the  available 
kinetic  data  to  the  pure  liquid  decomposition,  and  to  inquire  into  the 
relationship  between  the  low  temperature  kinetic  behavior  and  the  very 
fast  high  temperature  kinetics  associated  with  detonation. 

Extrapolation  of  kt  from  Eq.  (4-5)  to  the  rate  in  the  pure 
liquid  at  153°  indicates  that  t1j2  ^  200  hrs,  whereas  the  experimental 
value  (see  Section  B-l-a)  was  about  20  hrs.  However  this  latter  value  of 
tx/2  was  associated  with  the  period  of  rate  acceleration,  and  the  dis¬ 
crepancy  may  be  due  in  part  to  the  rapid  onset  of  catalysis  which  is  a 
consequence  of  the  higher  concentration. 

Since  the  decomposition  reaction  is  polar  and  shows  a  solvent 
effect,  perhaps  the  fundamental  question  concerns  how  nitrobenzene  and 
1,2  DP  differ  as  solvents  for  this  reaction.  At  this  time  we  can  do  no 
more  than  indicate  some  possible  differences  and  possible  consequences 
of  these  differences.  A  reaction  involving  a  polar  transition  state, 
i.e.,  a  state  having  charge  separation,  will  display  marked  sensitivity 
of  the  rate  to  the  dielectric  character  of  the  medium.21  Nitrobenzene 
has  a  high  dielectric  constant22  of  35.  The  value  for  1,2  DP  is  unknown; 
a  reasonably  close  model  may  be  1 , 2-dichloroethane ,  which  has  a  value22 
of  about  10.  On  this  basis,  one  would  expect  a  polar  reaction  to  proceed 
less  rapidly  in  1,2  DP  than  in  nitrobenzene,  contrary  to  the  observed 
result . 


If  hydrogen  fluoride  catalyzes  the  decomposition  by  interaction 
with  the  leaving  fluoride,23  then  the  ability  of  each  solvent  to  complex 


21A.  Streitweiser ,  Chem.  Rev.,  5(>,  615  (1956) 

22 J.  Hine,  Physical  Organic  Chemistry,  McGraw-Hill,  New  York,  1962,  p.  39 
23A.  Streitweiser,  Op.  Cit .  p.  622 
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with  HF  will  have  a  marked  kinetic  effect.  Nitrobenzene  can  probably 
complex  HF  effectively  through  interaction  of  the  aromatic  ring  with  the 
proton.  Similar  complexing  by  1,2  DP  would  be  expected  to  lead  to  de¬ 
composition  of  1,2  DP.  If  this  is  correct,  the  rate  of  decomposition 
should  be  faster  in  the  neat  liquid  than  in  nitrobenzene  solvent,  con¬ 
sistent  with  the  observed  results.  Unknown  factors  may  further  compli¬ 
cate  the  situation.  Until  more  is  known  regarding  the  mechanism  of 
catalysis,  we  cannot  properly  weigh  all  of  the  factors  affecting  the  change 
in  rate  in  going  from  solution  to  the  pure  liquid. 

Should  rate-controlling  steps,  alternative  to  the  one  proposed, 
become  predominant  at  high  temperatures  associated  with  the  detonation 
wave  (1000°-2000°K)  then  the  kinetic  parameters  derived  from  rate  studies 
at  176-187°C  become  inapplicable.  We  find  no  evidence  for  the  presence 
of  a  competing  unimolecular  or  bimolecular  reaction  based  on  product  dis¬ 
tribution  changes  on  changing  the  reaction  solvent  or  phase.  Nevertheless, 
a  reaction  which  proceeds  at  one-tenth  the  rate  of  the  principal  reaction 
at  187°C  (460°K)  will  constitute  26  percent  of  the  over-all  decomposition 
at  1000°K,  assuming  temperature-independent  frequency  factors  for  both 
reactions.  A  self-initiated  free-radical  chain  reaction  might  become 
important  at  high  temperatures,  despite  little  or  no  evidence  for  its 
presence  at  low  temperatures. 

Another  apparent  objection  to  the  use  of  low  temperature  kinetic 
data  for  detonation  phenomena  in  the  case  of  1,2  DP  arises  fran  the  fact 
that  the  final  product  of  low  temperature  decompositions  is  a  polymer, 
whereas  one  expects  gaseous  products  in  a  detonation  wave.  However  since 
most  organic  compounds  decompose  at  appreciable  rates  above  600°K, 
there  is  a  ceiling  temperature  above  which  no  polymer  is  formed;  thus, 
secondary  reactions  which  lead  to  polymer  at  low  temperature  are  supplanted 
at  high  temperature  by  alternative  decomposition  paths  which  give  gaseous 
products  and  probably  carbon. 
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C .  Kinetics  of  the  Thermal  Decomposition  of  2,2  DP 

1 .  Experimental 

Samples  of  2,2  DP  were  supplied  as  a  25  percent  solution  in  Aroclor. 
One  to  2  ml  quantities  were  distilled  on  the  vacuum  line  by  pumping  the 
vapor  through  a  nitrogen-cooled  trap  at  10~3  mm  Hg.  This  procedure  failed 
to  remove  all  of  the  2,2  DP  from  the  solvent,  even  though  we  have  found 
2,2  DP  to  be  rather  fugitive  in  nitrobenzene. 

Studies  of  the  decomposition  of  the  vapor  were  carried  out  in  the 
same  way  as  the  vapor  phase  studies  of  1,2  DP.  In  all  cases  neat  2,2  DP 
was  prepared  and  handled  remotely  on  the  vacuum  line.  Decomposition  at 
200°C  was  accompanied  by  formation  of  a  brown  polymer  on  the  walls  of 
glass  vessels;  some  evidence  of  etching  was  noted.  The  intermediate 
products  observed  in  the  gas  chromatography  were  different  from  those 
found  in  the  decomposition  of  1,2  DP.  A  typical  chromatogram  for  the 
decomposition  mixture  is  shown  in  Fig.  21. 

2 .  Results  and  Discussion 

The  few  vapor  phase  experiments  with  2,2  DP  at  200°C  indicate  that 
it  is  more  stable  than  1,2  DP  by  a  factor  of  five  or  more;  the  half-life 
of  2,2  DP  at  200°C  in  new  glass  is  about  57  hours.  Studies  in  nitro¬ 
benzene  solution  at  187°C  indicate  that  it  is  more  stable  in  solution  as 
well  . 


This  difference  in  rate  of  decomposition  must  reflect  a  difference 
in  mechanism,  with  initial  loss  of  NF2  or  HNF2  rate-controlling  in  the 
case  of  2,2  DP.  This  loss  of  HNF2  is  less  exothermic  than  loss  of  HF. 
Some  HF  was  formed  in  the  decomposition  of  2,2  DP.  One  bulb  heated  for 
80  hours  lost  0  35  x  10~3  moles  of  2,2  DP  and  gained  0.07  x  10"3  moles 
of  water.  This  corresponds,  assuming  that  all  of  the  HF  is  converted  to 
water  on  the  walls,  to  10  percent  of  the  HF  expected  if  each  molecule  of 
2,2  DP  loses  one  molecule  of  HF. 
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FIG.  21  CHROMATOGRAM  OF  DECOMPOSITION 
PRODUCTS  OF  NITROBENZENE 
SOLUTION  OF  2,2  DP  OBTAINED  ON  A 
5-FT  DINONYLPHTHALATE  COLUMN 
AT  76  C  AND  30  ml  min  FLOW  RATE 

Direct  formation  of  HF  with  concomitant  formation  of  N-f luoroimine 
is  improbable  for  steric  reasons.  A  more  plausible  path  for  HF  formation 

is 

4  HNF2  -♦  4HF  +  N2  +  N2F4 

with  a  possible  mechanism  involving  fluoronitrene  (NF): 


hnf2 

-> 

HF  + 

NF 

2NF 

N2Fa 

NF  + 

hnf2 

— > 

NF2  +  HNF 

2HNF 

— > 

Na  + 

2HF 

2NF2 

— ► 

NaF4 

2NaFa 

—¥ 

n2  +  N2F4 
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D .  Future  Plans 

Experiments  will  be  carried  out  to  determine  the  mechanisms  of 
catalysis  found  in  the  decomposition  of  both  1,2  DP  and  2,2  DP.  The 
structural  features  of  the  final  residues  from  both  these  reactions  will 
be  established  and  the  catalytic  activity  of  each  of  these  residues  will 
be  measured. 

A  more  rigorous  proof  will  be  provided  for  2- (f luoroimino)propioni tri le , 
one  of  the  decomposition  products  from  1,2  DP.  The  role  of  this  product 
in  the  over-all  decomposition  will  be  examined.  Similar  experiments  will 
be  made  for  intermediates  arising  from  the  d ecomposition  of  2,2  DP. 

Further  vapor-phase  decomposition  studies  will  be  carried  out  with 
both  1,2  DP  and  2,2  DP  at  higher  temperatures,  using  manometric  techniques, 
if  possible,  both  to  establish  the  kinetic  characteristics  of  these  re¬ 
actions  and  to  search  for  evidence  for  new,  high  temperature  reactions 
which  may  be  important  in  detonation  reactions. 

Studies  similar  to  those  described  in  this  report  will  be  done  with 
the  two  other  propane  isomers,  1,3  DP  and  1,1  DP  in  order  to  provide  as 
broad  as  possible  basis  for  generalizations  concerning  the  relation  of 
structure  and  reactivity  for  NF  compounds. 
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Appendix 

PREPARATION  AND  PROPERTIES  OF 
1,2-  AND  2 , 2-B IS(DIFLUOROAMINO)PROPANE 
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Appendix 

PREPARATION  AND  PROPERTIES  OF 
1,2-  AND  2 , 2-BIS (DI FLUOROAM I NO ) PROPANE 

Three  pounds  of  1 ,2-bis(difluoroamino)propane  were  prepared  by  the 
duPont  Eastern  Laboratories,  Gibbstown,  New  Jersey,  by  the  addition  of 
N2F4  to  propylene.  The  compound  was  shipped  as  a  20$  solution  in  methylene 
chloride.  Material  of  99$  purity  or  better  was  obtained  from  this  solu¬ 
tion  by  stripping  off  the  methylene  chloride  by  remote  distillation. 

Further  purification  on  a  dinonylphthalate  gas  chromatograph  column  pro¬ 
vided  material  of  99.8$  purity  for  the  kinetic  experiments. 

Approximately  1  pound  of  2 , 2-bis (dif luoroamino)propane  has  been  pre¬ 
pared  by  Aerojet-General  Corporation  in  Azusa  by  addition  of  difluoro- 
amine  to  acetone  in  sulfuric  acid.  The  compound  was  received  as  a  2  to 
1  solution  of  the  propane  in  Aroclor  1248.  Extensive  tests  of  the  ship¬ 
ping  solution  have  been  made  at  Aerojet  in  order  to  determine  if  any 
serious  shipping  hazards  existed.  Impact  sensitivity  values  for  the 
solution  were  reasonable  and  friction  tests  were  negative.  Static  sensi¬ 
tivity  tests  indicated  that  a  fire  hazard,  but  no  explosion  hazard,  existed 
for  the  2  to  1  solution.  However,  1  to  1  solutions  will  explode  by  spark 
initiation.  The  pure  material  is  sensitive  to  static  discharge,  being 
explosive  over  the  range  of  energy  that  a  human  body  can  accumulate. 
Isolation  of  the  pure  compound  (>99$)  was  obtained  by  remote  distillation 
under  reduced  pressure. 

Physical  properties  of  the  two  materials  are  summarized  in  Tables 
A-I  and  A-II. 
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Table  A-I 

PHYSICAL  PROPERTIES  OF 
1 , 2 -B I S ( D I FLUOROAM INO ) PROPANE 


Compound :  CH3  CH ( NF2 ) CH2 ( NF2 ) 
Boiling  Point:  B^C1 


Vapor  Pressure1 


T(°C) 

P(mm  Hg) 

-25  A 

4.2 

-13.0 

9.15 

0.1 

20.70 

10.4 

36.30 

24.4 

72.25 

Impact  Sensitivity 

(01  in  Mathieson  Machine) 


50  kg/cm1 
50  kg/cm2 
37  kg/cm2 
7  inches/5  kg4 


Static  Sensitivity4 

25,600  man  equiv,  volts 

Body  static  discharge  -  11,000  mev 


Freezing  Point:  -90°C 

Density:  1.256  g/ml  at  24°C;1  1.178  g/ml  at  25°C2 

Adiabatic  Compression  Sensitivity:  13  kg-cm/ml 

Nitromethane  10.4  kg-cm/ml1 

Thermal  Stability:  251°C  (decomposition  temperature);1 

Hot  Bar:  Evap'd.  at  250°C4 

Est’d  Activation  Energy:  73  kcal/mol1 

Heat  of  Formation:  -40.4  kcal/mole3 

Heat  Capacity:  0.348  ±  0.005  cal/deg/g5 


Stauffer  Chemical  Company,  Private  communication 
2Stanford  Research  Institute 

3Thermochemical  Panel  Meeting,  July  1961,  at  Aerojet-General 
Corporation,  Scaramento,  California 
4DuPont,  Private  communication 
SH.  Edlng,  Stanford  Research  Institute 
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Table  A-II 

PHYSICAL  PROPERTIES  OF 
2 , 2-BIS (DIFLUOROAMINO)PROPANE 

Compound :  CH3  C ( NF2 ) 2  CH3 

Boiling  Point:  73°^  (Vapor  press  at  25°  =  90  mm2) 

Freezing  Point:  About  0°^ 

Impact  Sensitivity:  53  cm/2  kg1 

"Erratic”3 

25-28  kg  cm  for  7  shots2 

Spark  Sensitivity:  Explosion  by  0.05  joule  discharge3 

(Body  discharge  energy — 0.02-0.07  joule)3 

Adiabatic  Compression  Sensitivity:  ’’Good”3 

Thermal  Stability:  Stable  "at  least  to  its  boiling  point"1 


1 Aero jet -Genera Y  Corporation ;  Report  0235-01-10 
aRohm  and  Haas;  Report  P60-24 

3J.  R.  Fisher,  Aerojet-General  Corporation,  Private  Communication 
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